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Abstract 

By using chaos expansion into multiple stochastic integrals, we make a wavelet analysis of two self-similar 
stochastic processes: the fractional Brownian motion and the Rosenblatt process. We study the asymptotic 
behavior of the statistic based on the wavelet coefficients of these processes. Basically, when applied to a 
non-Gaussian process (such as the Rosenblatt process) this statistic satisfies a non-central limit theorem even 
when we increase the number of vanishing moments of the wavelet function. We apply our limit theorems 
to construct estimators for the self-similarity index and we illustrate our results by simulations. 
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1. Introduction 

The self-similarity property for a stochastic process means that scaling of time is equivalent to an appro- 
priate scaling of space. That is, a process (l^t)t>o is self-similar of order > if for all c > the processes 
{Xct)t>o and {c^Yt)t>Q have the same finite dimensional distributions. This property is crucial in applica- 
tions such as network traffic analysis, mathematical finance, astrophysics, hydrology or image processing. 



We refer to the monographs j7[, [ij] or [28| for complete expositions on theoretical and practical aspects of 
self-similar stochastic processes. 

The most popular self-similar process is the fractional Brownian motion (fBm). Its practical applications 
are numerous. This process is defined as a centered Gaussian process {B^)t>o with covariance function 

Rnit, s) = E(i?f ) = i {t'" + s^" - \t - .p^) , t,s>Q. (1) 

It can be also defined as the only Gaussian self-similar process with stationary increments. Recently, this 
stochastic process has been widely studied from the stochastic calculus point of view as well as from the 
statistical analysis point of view. Various types of stochastic integrals with respect to it have been introduced 
and several types of stochastic differential equations driven by fBm have been considered. Other stochastic 
processes which are self-similar with stationary increments are the Hermite processes (see 12 1, an 



Hermite process of order q is actually an iterated integral of a deterministic function with q variables with 
respect to the standard Brownian motion. These processes appear as limits in the so-called nonccntral limit 
theorem and they have the same covariance as the fBm. The fBm is obtained for q = 1 and it is the only 
Gaussian Hermite process. For q = 2 the corresponding process is known as the Rosenblatt process. Although 
it received a less important attention than the fractional Brownian motion, this process is still of interest 
in practical applications because of its self-similarity, stationarity and long-range dependence of increments. 
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Actually the great popularity of the fractional Brownian motion in practice (hydrology, telecommunications) 
is due to these properties; one prefers fBm rather than higher order Hermite process because it is a Gaussian 
process and the corresponding stochastic analysis is much easier. But in concrete situations when empirical 
data attest to the presence of self-similarity and long memory without the Gaussian property (an example 
is mentioned in the paper jsol). one can use a Hermite process living in a higher chaos, in particular the 
Rosenblatt process. 

When studying self-similar processes, a question of major interest is to estimate their self-similarity order. 
This is important because the self-similarity order characterizes in some sense the process: for example in 
the fBm case as well as for Hermite processes this order gives the long-range dependence property of its 
increments and it characterizes the regularity of the trajectories. Several statistics, applied directly to the 
process or to its increments, have been introduced to address the problem of estimating the self-similarity 
index. Naturally, parametric statistics (exact or Whittle approached maximum likelihood) estimators have 
been considered. But to enlarge the method to a more general class of models (such as e.g. locally or 
asymptotically self-similar models), it could be interesting to apply semiparametric methods such as wavelets 
based, log-variogram or log-periodogram estimators. Information and details on these various approaches 
can be found in the books of Beran [7| and Doukhan et al. [13]. 

Our purpose is to develop a wavelet-based analysis of the fBm and the Rosenblatt process by using 
multiple Wiener-Ito integrals and to apply our results to estimate the self-similarity index. More precisely, 
let V' : R — >■ K be a continuous function with support included in the interval [0, 1] (called "mother wavelet"). 
Assume that there exists an integer Q > 1 such that 



tPi/'(i)dt = Oforp = 0, 1,...,Q-1 (2) 



and 



JtS. 



The integer Q > 1 is called the number of vanishing moments. For a stochastic process (Xt)tg[o 7v] and for 
a "scale" a G N* we define its wavelet coefficient by 



d{a,i) = ^ i' ( - - i] Xtdt = ^/a ijj{x)Xa(x+i)dx 
VaJ-oo \a. J Jo 

for i = 1, 2, . . . , Na with Na — [N/a] — 1. Let us also introduce the normalized wavelet coefficient 



(3) 



da,i = (4) 

(Ed2(a,i))5 

and the statistic 

VN{a) = —Y,[d^a,^)-^). (5) 

^ i—l 

The wavelet analysis consists in the study of the asymptotic behavior of the sequence Viv(a) when N oo. 
Usually, if X is a stationary long-memory process or a self-similar second-order process, then ¥.d'^{a,i) is a 
power- law function of a with exponent 2H — 1 (when a — > oo) or 2H + 1 respectively. Therefore, if Viv(a) 
converges to 0, a log- log-regression of -/^ X^tTi '^^('^i' *) onto aj will provide an estimator of H (with an 
appropriate choice of (aj)j). Hence, the asymptotic behavior of Vat (a) will completely give the asymptotic 
behavior of the estimator (see the Section [5] for details). There are four main advantages to use such an 
estimator based on wavelets: firstly, it is a semiparametric method that can be easily generalized. Secondly, 
it is based on the log-regression of the wavelet coefficient sample variances onto several scales and the graph of 
such a regression provides interesting information concerning the goodness-of-fit of the model (x^ goodness- 
of-fit have been defined and studied in Q or 0] from this log- regression) . Thirdly, it is a computationally 
efficient estimator (this is due to the Mallat's algorithm for computing the wavelet coefficients). Finally, it 



2 



is a very robust method: it is not sensitive to possible polynomial trends as soon as the number of vanishing 
moments Q is large enough. 

This method has been introduced by Flandrin [3| in the case of fBm and later it has been extended to 
more general processes in [l| . The asymptotic behavior of the estimator obtained from the wavelet variation 



has been specified in the case of long- memory Gaussian processes in [6( or |21| . for long memory linear 



processes in [27| and for multiscale fractional Gaussian processes in . However the case of the Rosenblatt 
process has not yet been studied (note that the wavelet synthesis of Rosenblatt processes was treated in Q 
and we will use this method in the section devoted to simulations). 

In our paper we will use a recently developed theory based on Malliavin calculus and Wiener-Ito multiple 
stochastic integrals. Let us briefly recall these new results. In ^] the authors gave necessary and sufficient 
conditions for a sequence of random variables in a fixed Wiener chaos (that means in essence that these 
random variables are iterated integrals of a fixed order with respect to a given Brownian motion) to converge 
to a standard Gaussian random variable (one of these conditions is that the sequence of the fourth order 
moments converges to 3 which represents the moment of order 4 of a standard Gaussian random variable). 



Another equivalent condition is given in the paper 2J] in terms of the Malliavin derivative. These results 



created a powerful link between the Malliavin calculus and limit theorems and they have already been used 
in several papers (for example in 33 1, and 10 1 to study the variations of the Hermite processes). 



Recall (see Q and 13) that if X = is a fBm in ([5|) then the following fact happens: for any Q > 1 
and H G (0, 1) the statistic VAr(a) renormalized by VlV converges to a centered Gaussian random variable. 
If Q = 1 then the barrier H — i/A appears: that is, the behavior of yjv(a) is normal {i.e. it satisfies a central 
limit theorem) only if G (0, 3/4) and we prove in Section [3] that the limit of VAr(a) (normalized by iV^^^^) 
is a Rosenblatt random variable when H E (3/4, 1). This non-central limit theorem is new in the case 
of the wavelet based statistic (it is known for other type of statistics, such as the statistic constructed from 



the variations of the process, see 33|). In this case we also prove that the limit holds only in law and not in 
in contrast to the case of quadratic variations studied in [s^ • 



The study of Vat (a) in the Rosenblatt process case (see Section |3l formula ([TT]) for the definition of the 
Rosenblatt process) with H G (1/2, 1) put in light interesting and somehow intriguing phenomena. The main 
fact is that the number of vanishing moments Q does not affect its convergence and the limit of Vn (a) is 
always non-Gaussian (it is still Rosenblatt). This fact is unexpected and different to the situations known in 
the literature. Actually, the statistic V^v can be decomposed into two parts: a term in the fourth chaos (an 
iterated integral of order 4 with respect to a Wiener process) and a term in the second chaos. We analyze 
here both terms and we deduce that the term in the fourth Wiener chaos keeps some of the characteristics 
of the Gaussian case (it has to be renormalized by Vn and it has a Gaussian limit for H G (1/2, 3/4)). But 
the main term in V/v(a) which gives the normalization is the second chaos term and its detailed analysis 
shows that the normalization depends on H (it is of order of N^^^) and its limit is (in law) a Rosenblatt 
random variable. 



The consequences of these results are also interesting for the wavelet based estimator of the self-similarity 
order. Assume that a sample {Xi, ■ ■ ■ ,Xn) is obseved, where X is a fBm or a Rosenblatt process. First, 
by approximating the wavelet coefficients ([3]) , we consider a statistic V^v (a) that can be computed from the 
observations {Xi, ■ ■ ■ ,Xn) and we prove that the limit theorems satisfied by Vjv(a) also hold for VN{a) as 
soon as a is large enough with respect to N. Secondly, we deduce the convergence rates for the wavelet based 
estimator of H following the cases: Q > 2 and X is a fBm, Q = 1, H G (3/4, 1) and X is a fBm or H E (i, 1) 
and A" is a Rosenblatt process. The regularity of also plays a role. For practical use, it is clear that if 
A is a fBm it is better to chose Q > 2 and the mother wavelet i/j to be twice continuously differentiable 
(for example, this is the case when is a Daubechies wavelet with order > 8). On the other hand, if X 
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is a Rosenblatt process the number of vanishing moments Q plays no role. Our simulations illustrate the 
convergence of ¥^{0.) and of the estimator of H in this last case. 

Our results open other related questions: for a stationary long-memory Gaussian or linear process, para- 
metric estimators (such as the Whittle's estimator, see [l?]] and [18|) or semiparametric estimators (such as 



wavelet based estimator, see 



21 



and 



27| ) satisfy a central limit theorem with an "usual" convergence rate 



(that is \/N for Whittle's estimator and yjN/a for wavelet based estimator). This is not the case of the 
Whittle's estimator for non-linear functionals of a long-memory Gaussian process (see 19|). As we prove 
in our paper, this is no longer the case for the wavelet estimator based on the observation of a Rosenblatt 
process. Therefore it seems that the Rosenblatt processes and the second order polynomials of Gaussian 
process give the same asymptotic behavior for the estimators. Indeed, in 19| it has been established that the 
Whittle's estimator in the case of second order polynomials of Gaussian processes satisfies a noncentral limit 
theorem with a Rosenblatt distribution as a limit and with the same convergence rate as in Theorem|3]below. 
Several interesting questions arise naturally from this: firstly, are there any reciprocal results true? (i.e. may 
we suspect that the Whittle's estimator of the long-memory parameter for the increments of the Rosenblatt 
process and the wavelet based estimator of the long-memory parameter for a second order polynomial of 
a long-memory Gaussian process also satisfy the same noncentral limit theorem?). Secondly, are the con- 
clusions the same for other semiparametric estimators such as log-periodogram or local Whittle's estimators? 



We organized the paper as follows. Section [2] contains some preliminaries on multiple Wiener-Ito integrals 
with respect to the Brownian motion. In Section[3]wc treat the situation when the driven process is the fBm. 
In this case our new result is the noncentral limit theorem satisfied by the wavelet based statistic proved 
in Theorem 2. In Section |4] we enter into a non-Gaussian world: our observed process is the Rosenblatt 
process and using the techniques of the Malliavin calculus and recent interesting results for the convergence 
of sequence of multiple stochastic integrals, we study in details the sequence V/v(a). In Section[5]we construct 
an observable estimator based on the approximated wavelet coefficients and we study its asymptotic behavior. 
We also construct an estimator of the self-similarity order and we illustrate its convergence by numerical 
results. 



2. Preliminaries 

2.1. Basic tools on multiple Wiener-Ito integrals 

Let (VFt)te[o,T] be a classical Wiener process on a standard Wiener space {^l,F,P). If / S L^([0, T]") 
with n > 1 integer, we introduce the multiple Wiener-Ito integral of / with respect to W. We refer to [23 1 
for a detailed exposition of the construction and the properties of multiple Wiener-Ito integrals. 



Let / G 5„, that means that there exists n > 1 integer such that 

^ ^ Cii,...,i„lAi. x...xAi„ 

where the coefficients satisfy Ci^^...^i^ = if two indices ik and ii are equal and the sets Ai E B{[0,T]) are 
disjoint. For such a step function / we define 

Uf)^ E c,,,...,,„w{A,)...w{AJ 

where we set W{[a,b]) — Wb — Wa. It can be seen that the application /„ constructed above from 5„ 
equipped with the scaled norm -^ \\ ■ ||l2([o.t]") to L^(f7) is an isometry on 5„, i.e. for m, n positive integers, 

E(/„(/)/m(g)) = '^!(/).9>L2([o,T]") ifm^n, (6) 
E(/„(/)/,„(5)) - ifm^n. 
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It also holds that 

Uf)^In{f) (7) 

where / denotes the symmetrization of / defined by /(xi, . . . ,a;„) = ^ J2a£S„ fi^<y{i)^- ■ • ,a;cr(n))- 
Since the set 5„ is dense in L^([0, T]") for every n > 1, the mapping /„ can be extended to an isometry from 
L^([0,T]") to L^(r2) and the above properties ^ and ([7]) hold true for this extension. Note also that /„ can 
be viewed as an iterated stochastic integral 

Inif) =n\ I r ■■ I' • • ■^in)dWt, . ..dWt^. 



Jo Jo 



We recall the product formula for two multiple integrals (see [23|]): if / e L^([0,r]") and g e L^([0,r]™) 
are symmetric functions, then it holds that 



m/\n 



Unimig) = E eiCi,Cil^+n-2iif ®i g) (8) 

t=0 

where = "l^^i and the contraction f ®ig belongs to L^([0, 7"]™+"-2^) for £ = 0, 1, . . . , to A n with 

(/ ®i • ■ • , Sn-t,ti, . . . ,tm-e) = / /(si, . . . , s„-e, ui, . . . ,ut)g{ti, . . . ,tm-e,ui, ui)dui . . . dui. 

J[q,tY 

In general the contraction f <S)£ g is not necessarily a symmetric function even if the two functions / and g 
are symmetric. 

3. The case of the fractional Brownian motion 

3.1. A presentation using chaos expansion 

We will assume in this part that X = a (normalized) fractional Brownian motion (fBm in the sequel) 
with Hurst parameter H G (0, 1). Recall that is a centered Gaussian process with covariance function 
([1]). It is the only normalized Gaussian i7-sclf-similar process with stationary increments. Recall also the 
fBm (i3^)tg[o.Ar] with Hurst parameter H e (0, 1) can be written 



= / K"{t,s)dWs, te[0,N] 
Jo 



where {Wt,t € [0,A^]) is a standard Wiener process and for s < t, and for _ff > i the kernel has the 
expression 

K"{t, s) = CHS^^" / ~ s)"~^u"-idu (9) 

with ch = ( ^(2^^2g g- j-) ) ^ ^^'^ beta function. For t > s and for every H e (0, 1) the partial 

derivative of with respect to its first variable is given by 

{t, s) = diK"{t, s) = CH ('^Y'" [t - s)^^i (10) 



dt ' ' ' ' ' ' ' " \ t 

In this case it is trivial to decompose in chaos the wavelet coefficient d{a, i). By a stochastic Fubini theorem 
we can write 

d(a,z) = ^/^ j^^{x)B^^.^^i^dx^^ j^4,{x)dx(^j^^ ^ dB^^ 

= f ^{x)dx f ^ '^ \"ia{x + i),u)dWu^Iiifa4-)) 

Jo Jo 
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where Ii denotes the multiple integral of order one (actually, it is the Wiener integral with respect to W) 
and ^ 

/a,i(M) = l[o,Q(i+i)](M)\/a / ilj{x)K"{a{x + i),u)dx. (11) 

J(f -i)VO 

Thus, for all a > and i e N, 

with G4,{H) = -l [ [ ^p{x)^{x')\x- x'f"dxdx' (12) 
2 Jo Jo 

(see [1]). This formula is essential for the estimation of H (see Section [5|). Using the product formula ([8]) 



h{f)h{g) ^h{f®g) + (/,ff)L2([o,w]) 

we obtain 

(Ed(a,z))2 lj-^2(/^ ) 

where 

/(;)^a-^-iC,(i/)-i-^^/,«,^. (13) 

I— 1 

As a consequence the wavelet statistic Vn can be expressed as a multiple Wiener-Ito integral of order two 
and therefore the chaos expansion techniques can be applied to its analysis. 



3.2. A multidimensional Central Limit Theorem satisfied by {VN{ai))^ 



<i<e' 



When the observed process is a fBm with H < 3/4, the statistic V/v(a) satisfies a central limit theorem. 
This fact is known and we will not insist on this case. We just recall it to situate it in our context. Since 
^^iif) = 2!||/||l2[o Ar]®2, if {ai)i<i<i is a family of integer numbers such that ai = ia iov i — I, . . . ,£ and 
a e N*, it holds that 



1 1 



Cov(VW(a,),Viv(a,)) = 2\{pqa'r^"''C^{H)-^ — C-f^^^iioM)^ 



j=l j' = l 

We know (see e.g. Q and Q) that: 

(/Qpj,/a,j')L2([o,Ar]) = ^{d[ap,i)d{aq,i')) 

= -]-{pqa^Y'^a^" [ [ '4,{x)i}{x')\px^qx' +pj-qj'\^"dxdx'. (14) 
2 Jo Jo 

and using the Taylor expansion and the property ^ satisfied by if), 

=^ ICov (yA.(ap), VM{a,))\ < C ^ 0(l + - gj'!)'^' 



r-l /-l 



iV2 

i=l j' = l 



Consequently, if Q > 1 and G (0, 1) or if Q = 1 and H e (0, 3/4), 
Na Gov {VN{ap),VN{ag)) — ^ h{p,q,H) with 

4(., ,.H) = j^^^j^ f f ^(.)^(.')bx - + kdp,rdxdx')\ (15) 

where dp^ is a constant depending only on p and g. Moreover, the following limit theorem holds true. 
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Theorem 1. Let {Xt)t>o be a fBm, Vn^u) be defined by and Lf \H) — {£i{p,q, H))Kp q<i with 
£i{p, q, H) defined in U^l . Then if Q > 1 and H G (0, 1) or if Q ^ 1 and H G (0, 3/4), for all a > 0, 

^/KiVN a)) Mm (0 , if ^ {H)) . (16) 

Proof: It is well-known in the literature (see e.g. [sj). ■ 

The above central limit theorem does not hold for Q — 1 and iJ G (|, !)■ We will focus on this case in the 
following paragraph. 

3.3. A noncentral limit theorem satisfied by yjv(a). 

We need at this point to define the Rosenblatt process. The Rosenblatt process with time interval [0, N], 
denoted in the sequel by {Rt^)te[o,N] appears as a limit in the so-called noncentral limit theorem (see [12| 



or 



29[). It is not a Gaussian process and can be defined through its representation as a double iterated 



integral with respect to a standard Wiener process (see [32] )■ More precisely, the Rosenblatt process with 
self-similarity order He (5,1) is defined by 

^ f f Lf{y,,y2)dWy,dWy, (17) 
Ja Jo 

where {Wt, t G [0, A'^]) is a standard Brownian motion and 

Xf(yi,y2) -d//l[o,t] (2/1)1 [o,t] (2/2) / diK"'{u,yi)diK'''{u,y2)du, (18) 

J aiVj/2 

with the standard kernel defined in (jH]) appearing in the Wiener integral representation of the fBm, its 
derivatives being defined in (jlOp and 

, H + 1 1 / H , , 

The random variable R^ is called a Rosenblatt random variable with self-similarity index H. A Rosenblatt 
process is a process having stationary increments and 

• it is iJ-self-similar in the sense that for any c > 0, (i?^)f>o [c^ R^)t>o, where " " means 

equivalence of all finite dimensional distributions; 

• '&{\R^\'P) < 00 for any p > 0, and R^ has the same variance and covariance as a standard fractional 

Brownian motion with parameter H. 

• the paths of the Rosenblatt process are Holder continuous of order 5 < H . 

We obtain the following noncentral limit theorem for the wavelet coefficient of the fBm with H > j. Define 



^^("^^ - I 4ff-3 ) C^) =d^, CW) ■ ^ ^ 

Then, 

Theorem 2. Let {Xt)t>o be a fBm and Viv(a) be defined by Lf Q = 1 and j < H < 1, 



where LIq = 2H — 1, £2{H) is defined by iSOp and R^° is a Rosenblatt random variable given by ( |_?7| ). 
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Proof: With /^"' defined as in (jl3p . we can write 

Using the expression PT|) of the kernel fa^i, one has 
1 1 ^" 

/tv" (yi:2/2) = q2g(g^(jj-) ]^ Xll[n.«(''+i)l(2/i)l[o,a(j+i)](y2) 

X /" /" ip{x)ip{z)K" {a{x + i),yi)K" {a{z + i),y2)dxdz. 

7(!£L_i)vO J(^-j)VO 

To show that the sequence i2^{H)N^^^^ I2 (/^'') converges in law to the Rosenblatt random variable R^^~^ 
it suffices to show that its cumulants converge to the cumulants of (see e.g. [l3| where it has been 

proven that the law of a multiple integral of order 2 is given by its cumulants). On the other hand, we know 
(see 17 1, [2^) that the fc-cumulant of a random variable /2(/) in the second Wiener chaos can be computed 
as follows 

ckihif))^ / dyi...dykfiyi,y2)f{y2,y3)---fiyk-i,yk)f{yk,yi)- (21) 

"'[0,1]'= 

In particular the fcth cumulant of the Rosenblatt random variable R\^^^ is given by 

Ck{Rl"~^) = d^H-i{H{2H - 1))^ / [\xi - X2I • . . . • \xk-i -Xk\-\xk- xi\]^"'^dxi . . . dxk- 

J[o,i]fc 

We compute 

c.(c^(ff)7vr^^/2(/^"^)) 

= Ar_(2^f-2)fc^-fc / ^y^ _ _ / dxidzi . . . dxkdzk^l^{xi)'il;{zi)il;{x2)'ip{z2) ■ . . i;{xk)^{zk) 

X K"{a{xi + ii), yi)K"{a{zi + ii),y2)K" {a{x2 + ^2), y2)K" {a{z2 + ^2), ys) x . . . 

... X K"{a{xk-i +ik-i),yk-i)K"{a{zk + ik),yk)K" {a{xk + ik),yk)K" {a{zk +ik),yi)- 

Using Fubini Theorem and the fact that (for every x, x' , i, j, a) 

j-a{x+i) /\a{x' + j) 

/ K"{a{x + i),yi)K"{a{x' + j),yi)dyi = RH{a{x + i),a{x' + j)) 

Jo 

= a^" Rh{x + i,x' + j) 
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(from the representation of the fBm as a Wiener integral with respect to the Wiener process) we get 

= N^^^'^^'^N-'' V / dx^dz^...dxkdzkxl^{xi)^{zM{x2)^{z2)...xk{xkmzk) 
n,......=W[0-il^'= 

Rh{zi + ii,X2 + i2)RH{z2 + i2,X3 + is) ■ ..RHizk-i + ik-i,Xk + ik)RH{zk + ik.xi + ii) 
= iVi^^-'^^a^iffc^-fe / dxidzi...dxkdzktP{xi)ij{zi)tP{x2Mz2)...i^{xkMzk) 



'[0,1]^'' 

X [\zi - X2 + ii - i2 \ ■ \z2 - X3 + i2 - ial . ■ ■ ■ \zk-i - Xk + ik-i - ik\ ■ \zk - xi + ik - ii\ 



2H 



^2H 



[0,1]^'= 



ii,...,ifc = l 

dxidzi . . . dxkdzk^{xi)ip{zi)^{x2)tp{z2) ■ ■ ■ ■)p{xk)tp{zk) 



n - 12 



2H 



1 + 

ik - ii 



2H 



N(^"''^'^H\2H-\fN-'^ E {W-^2\■.■■\^k-l-^k\■\^k-^l\) 

ii,...,ik=l 



2H-2 



[0,1] 



dxidzi . . . dxkdzkip{xi)'ip{zi)Tl>{x2)ij{z2) . . . ^p{xk)ip{zk)xiZi . . . x^zu 



and we used the fact that the integral of the mother wavelet vanishes and a Taylor expansion of second order 
of the function (1 +x)^^. By a„ ^ fe„ we mean that the sequences a„ and 5„ have the same limit as n — > oo. 
As a consequence, by a standard Riemann sum argument, the sequence 



Na 



j^i2H-2)kj^-kJ2 (|,,_.,| 
ii,...,ifc = l 

converges to the integral 



X . . . X - lk\\lk ~ « 



ii,...,ik = l 



{\xi - X2I X . . . X \xk-i - Xk\\xk - xi\y^ ^dxi . . .dx 



|^l - Z2| X . . . X \ik-i - ^fc||^fc - ^l| 

Na 



, 2H-2 



[0,1]'= 



and therefore it is clear that the cumulant of £2 ^{H)N^ l2{f'N^) converges to the k cumulant of the 
Rosenblatt random variable RI"~^ (see 0, lil)- ■ 



In the case of the statistic based on the variations of the fBm, in the case H G (3/4,1) the statistic 
— E] '^j\j~2H " renormalized by a constant times iV^^^^, converges in L^(i7) to a Rosenblatt 

random variable at time 1 (see [s^). In the wavelet world, our above result gives only the convergence in 
law. The following question is then natural: can we get L^-convergence for the renormalized statistic V;v(a)? 
The answer is negative and a proof of this fact can be found in the extended version of our paper available 
on arXiv. But we will present below a brief and heuristic argument to see what the convergence does 
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not hold. The term f'^^ can be written as 

1 1 ^" 

In iyi,y2) = ^2ff^r^ ]y- ^ l[0,aCi+l)](2/l)l[0,a(»+l)](2/2) 



1 [0,ai] 

^ Ja Jo 

+ l[o,ai](yi)l[aj,a(i+j)](2/2) J j dxdzip{x)ijj{z)K" {a{x + i) ,yi)K" {a{z + i) , y2) 

+ l[o,ai](y2)l[a»,a(i+i)l(yi) j j dxdzijj{x)ijj{z)K" {a{x + i) ,yi)K" {a{z + i) , y2) 

+l[ai,a(i+»)](yi)l[M,a(i+»)](2/2) / / dxdzil){x)'4){z)K" (a(x + i),yi)K" {a{z + i),y2)^ 



fN^^\yi,y2) + fN''^\yi,y2) + fN'^\yi,y2) + .fN'^\yi,y2)- 



It can be shown that the terms Nl"^" f\^-'^\ N^'^" fl^'^^^ and N^-^" f*^'^^ converge to zero in h^{[0,oo)'^) 
as Na — >■ oo. It remains to understand the convergence of the term Z^'^"*. Using again the property 1^ of 
the mother wavelet tp we can write 

fN'^\yi,y2) = 2Hr^ (u\ i:rYl / / ^2;dzV'(x)V'(z)i[o,„](yi)i[o,«](y2) 

a ) i\a Jo Jo 

X {K" (aix + i),yi) - K" (ai, yi)) [K" {a{z + ^), 7/2) - K" (ai, ^2)) - 

Therefore, with a(a, i, a;) and /3(a, i, z) located in [ai, ax + ai] and [ai, az + ai] respectively, 

1 1 , /•! /•! 

hifN'^^) = (Tj\ V / / da;dzV'(x)V'(2;)l[o,„](2/i)l[o.«](y2) 

xaxdiK^ {a{a,i,x),yi) x az diK^ {I3{a,i, z),y2) 

By approximating the points the points Q;(a, i, x) and /3(a, i, z) by and since (from an usual approximation 
of a sum by a Riemann integral) when A^^ — > 00, with 2/1,1/2 G [0, A^], 

X! l[o,«](yi)l[o,a»](2/2)9i-fir^(aj,yi)9iii:^(ai,?/2) ^ / du diK^ {au,yi)diK" {au,y2) 

-■_ r> J ftyi \/V'>) /a 



i=0 



(j/lVi/2)/a 

02^-1 a 

where L'^^^^ is the kernel of the Rosenblatt process with self-similarity index 2H — 1 (see its definition 
in (US])), the sequence £2'^ iH)N^^'^^ f'j^'^'^ is equivalent (in the sense that it has the same limit point- 
wise) to N^-'^"L^jf-\ Then, in some sense, ^^^(i7)iV2-2ffy^(-^) equivalent to N'^'^" hih]^^^) = 
]\l^-2H j^2H-i ^ ^Yiis equivalence is only in law. The fact that the sequence N^~'^"VNia) is not 

Cauchy in comes from the fact that the sequence N^~'^^ R^~^ is not Cauchy in as it can be easily 
seen by computing the square mean of the difference N^^^^ R^^^ — M^^'^^ R^^^ . 



4. The Rosenblatt case 



4-.1. Chaotic expansion of the wavelet variation 

We study in this section the limit of the wavelet-based statistic Vn given by ([5|) in the situation when the 
observed process X is a Rosenblatt process. Throughout this section, assume that X = is a Rosenblatt 
process with self-similarity order H > 1/2 defined via the stochastic integral representation pT|) . We first 
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express the statistic (a) in terms of multiple stochastic integrals. The wavelet coefficient of the Rosenblatt 
process can be written as 



d{a,i) ~ \ 
where, with H' and cIh defined in ([T9l) . 



a(x-\'i) pa{x-\-i) 



La{x+^){yi.y2)dWy,dWy^ 



ga,i{jjl,y2) = dH ^/al\Q^a{i+l)]{yl)\o,a(i+l)]{y■. 



X I dxip{x) ( 

/(ii^_i)vo ^ 



a(x-\-i) 



yiVy2 



diK"\u,yi)diK"'{u,y2)du 



for every j/i,?/2 > 0. The product formula for multiple stochastic integrals (|8]) gives 

W)h{9) = W ®9) + 4/2(/ ®i g) + 2(/, (7)l2(o_^p2 
if f,g £ L^([0, A^]^) are two symmetric functions and the contraction f (E)i g is defined by 

U ®i g){yi^y2) = / fiyi,x)g{y2,x)dx. 
Jo 



Thus, we obtain 



d\a,l) = /4 (<7®?) +4/2(5a,»«)l <?a.0 + 2||5a,»|| 



2 



and note that, since the covariance of the Rosenblatt process is the same as the covariance of the fractional 
Brownian motion, we will also have 

Therefore, we obtain the following decomposition for the statistic VN{a): 



VNia) 



No. 



with 



T2 = a 
14 — a 



c„{H)-'±Y^=.h{gf}) 



(22) 



To understand the limit of the sequence Vn we need to regard the two terms above (note that similar terms 
appear in the decomposition of the variation statistic of the Rosenblatt process, see 33|). In essence, the 
following will happen: the term T4 which lives in the fourth Wiener chaos keeps some characteristics of 
the fBm case (since it has to be renormalized by \/Na except in the case Q = 1 and H > j where the 
normalization is N^^'^^) and its limit will be Gaussian (except for Q = 1 and H > |). Unfortunately, this 
somehow nice behavior does not affect the limit of Vn which is the same as the limit of T2 and therefore it 
is non- normal (see below). 

Now, let us study the asymptotic behavior of ET4 . From (|22|) . we have 



T.^hig^^^) where g^^^ = a-'"-'C^iH)-' 9^1 



Na 



(23) 
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and thus, by the isometry of multiple stochastic integrals, 



IT^ — 4!||.g^^|jL2[o^Ar]4 



where by g^j^^ we denoted the symmetrization of the function g]^' is its four variables. Since |157v^IIl2[o jv]* — 



l;;(o)||2 



llfl'VllL2[0,Af]4 



we obtain 



1 



But, 



< AIC^{H) ^— ^{ga,z,ga,j)l^Q^N]<S2 



{9aS,ga,j)i?(o,N]»^ = ^ E(d(a, i)d(a, j ) ) . 



Again, using the fact the fBm and the Rosenblatt process have the same covariance, we obtain the same 
behavior (up to a multiplicative constant) as in the case of the fractional Brownian motion. That is, using 
([TS]) and the proof of Theorem [2] 

If Q > I and 
(24) 



Proposition 1. Let {Xt)t = {Rf)t>o be a Rosenblatt process and T4 be defined by 
H e {\, 1) orifQ^l and H G (i, |), then with 1, H) defined in l[T5\) 



iV,E(r|) <7Va4!||5j;>||L2[o,w]4 ™> 3£i(l,l,i/) 



and, if Q — 1 and H € (|, 1) then with i2{H) defined in 



N- 



(25) 



The above result gives only an upper bound for the LP' norm of the term T4; this will be sufficient to obtain 
the desired limit of the sequence Vat (a) in the Rosenblatt case. It will follow from the following paragraph 
that the L^ norm of the term T4 will be dominated by the L^ norm of the term T2 and therefore the limit 
of T2 will be the limit of the statistics Vat (a). 

4.2. Asymptotic behavior of the term T2 

We study here the term in second Wiener chaos that appears in the decomposition of Vn (a) ■ We have 

(26) 



T2^l2{h^N^) with /ii;^=4. 



1 



We first compute the contraction ga^i (8)1 ga,i- We have 
{9a,t®iga,i){yi,y2) = / gaAyi^z)ga,tiy2,z)dz 



a dji l[o,a(i+i)] (2/i)l[o,a(i+i)] (2/2 



Q(4+1) 



a(x+z) 



yiVz 



dx 



a d\ l[o,a(j+i)] (2/i)l[o,a(j+i)] (2/2) ( 

fa{x-\'i) i*a{x' +1) 



Xx')( 



dz I dxtp{x 

diK"'{u', y2)diK"' {u' , z)du'^ 



diK"' {u,yi)diK" {u, z)du 



H' I 



y2Vz 



dx ip{x) 



(ia_i)vo 



dx' ^{x') 



VI 



AI{u,yi,u' ,y2) du du' / M{u, z,u' , z)dz 



V2 
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where M{u,yi,u' ,y2) — diK^ {u,yi)diK^ ("',2/2) and H' = {H + l)/2. Now, we have already seen that 
J*'^'^ {t, z)K^ {s, z)dz — Rfi{t,s) with Rfi{t^s) given in ([IJ and therefore (see [2^, Chapter 5) 

/ M{u,z,u',z)dz = H'{2H' ~\)\u-u'\'^"'~^ (27) 
Jo 

(In fact, this relation can be easily derived from /J"^" {u, yi)K^ [v, yi)dyi — Rh' {u, v), and will be used 
repeatedly in the sequel). Thus denoting aH — H'{2H' — 1) = H{H + l)/2 and since ^ is [0, l]-supported, 
we obtain 



{9a.t<E)l ga.i){yi,y2) = arfif^H l[0,a(j+l)](yi)l[0,a(j+l)](y2) / / dxdx' tp{x)ll;{x') 



u-u'\ M{u,yi,u',y2)dudu'. 

By direct computation, it is possible to evaluate the expectation of T| as follows (the proof can be found on 
the extended version on arXiv) 

iVr^^ETl ^ 32—-pi^-^(f ij{xMx')xx'\ux-vx'\''''-'dxdx'dudvy 



[0,1]" 



N^oo ' H{H +1)^ \C^{H) J V C^{H) J ' 



where we used / xx \ux — vx \ dudv— 



2H'{2H' - 1) 



\x\— + \x 



x'\'^^ ) and a\j d% = j {2H — l)^i7^. We do not prove here this estimate because it is a consequence of 
the following proposition which show that the sequence C^^{H)N^^^T2 (and therefore the sequence ¥^{0)) 
converges in distribution to a Rosenblatt random variable with self-similarity index H . 

Proposition 2. Let (i?^)f>o be a Rosenblatt process and let T2 be the sequence given by i2^) and computed 
from (i?^)t>o. Then, for any Q > 1 and H G {^^1), with Ct^ given by \28]). there exists a Rosenblatt 
random variable R^ with self-similarity order H such as 

C^l{H)Nl-"T2 A i?f. 

Proof: This proof follows the lines of the proof of Theorem [51 With T2 = l2(h!f^^) in mind, as in the 
proof of Theorem [21 a direct proof that cumulants of the sequence N^~^ l2[h^^^) converge to those of the 
Rosenblatt process can be given. Indeed, since the random variable iV^~^/2 (/ij^'') is an element of the 
second Wiener chaos, its cumulants can be computed by using formula (j2ip . By using the key formula (j27p 

c,(ivr^/2(4"^)) 

^N^^'-"^ [ dyi...dykh^;\yi,y2)h^^\y2,y3)x---xh''^\yk,yi) 
J[o,i]fc 

A', 



-Hk^-2Hk(j^jjyk^k^2k^2k f TT ^(^ .)^(:^;)d:,^.d:,; 

ra{xi+ii) ^a(x\+ii) ^a{xj+ij) ^a(x'-+ij) k s2_f/'-2/jL 

/ / . . . / / duidu'i . . . dujduj {\\\uj - "i+i I 

Jo Jo Jo Jo ^ .J ^ ^ 
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with the convention Ufc+i = Ui. Next, we will make the change of variables uj = auj and u'j = au'j and this 
will simplify the factors containing a. We thus get 

pxi+h i-x[+ii j-Xk+tk rx'i^+ik , J_ s2_f/'-2/jL. ^2H' -2 

X / / •••/ / duidu[...dukdu'ki^[[\uj -u'jlj (lim-Wj+ilj 

•'0 ^ ''0 i — 1 i — 1 



We will note at this point that the vanishing moment property of the wavelet function ^ allows us to replace 
integration intervals [0, xj + ij] by intervals [ij,Xj + ij]. After doing this, we perform the change of variables 
Uj = Uj — ij, u'j = u'j — ij (for every j = 1, . . . ,k) to obtain 

= N-"^C^{H)-HH^>}a'^ / Y{^j{xj)i,{x'j)dxjdx'j 

nx[ rxk rx'k k ^2ff'-2/jL n2H'-2 

"Jo Jo '^"i'^"'i---^"'='^"Hlll"j~"^lj \\{Wj-u'j+i+ij-ij+i\j 



X 



and then, with changes of variables Uj = ^ and u'j = we can write 

Na „ k 

= N-"''C^{H)-''A''dj^aj^ y / Y\xjx'.iP{xj)iP{x'^)dxjdx'j 



N-"''C^{H)-H'^dfal' Y / llxjx'ji^{xjMx',)dxjdx'j 



An analysis of the function (l+x)"^^ ~^ in the vicinity of the origin shows that the cumulant Ck{N^~^ ^{h'"^^)) 
behaves as (we recall that by a„ ~ &„ we mean that the sequences a„ and 6„ have the same limit as n ^ oo) 

CkiK-^hih'-^'^)) ~ N-"''C^iH)-H>^dfaf [ n i'{xjmx',)dx,dx'^{xjx'jr 

•^[0,1]"'= j = l 

k p k 2 if' 2 

ii,...,ifc = lj=l i[0,l]2fe j^-^ 

By using an usual Riemann sum convergence, we notice that 

Na k Na k . I s 2/^-2 

ii,...,ife = l j=l ii,...,ife = l j=l ^ o / 

converges as A^o — >■ oo to the integral 

/ . . . da;fc(|xi - a;2| • |a::2 - aJsl X • • • X - ^ 

y[o,i]fe 
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which is dj^aj^Ck{Ri) (here Ck{Ri) denotes the cumulant of the random variable Ri). We conclude that 



,2H'-2 



X / duidu'i . . .dukdu'i^l I I \ujXj — u'jx'A 

C^{H)-HH^H(^1j + {x'f"' -\X- X'\''"' dxdx'^ Cfc(i?f ) 



C4H)-H''d1jC^{H'fck{R?) 

where we recall the notation C^{H) — — ^ jj^ ~ x'p^ dxdx' . As a consequence C^^T2 con- 

c^(g') 



verges in law to R^ where Ct^ = Adn ^'*''^ ^ 



We finally state our main result on the convergence of the wavelet statistic constructed from a Rosenblatt 
process. Its proof is a consequence of Proposition [1] and [2] 

Theorem 3. Let {Xt)t>o = {R") t>o be a Rosenblatt process and Vat (a) be defined by Then, for any 

Q > 1 and G 1), there exists a Rosenblatt random variable R^ with self-similarity order H such as 

C^^{H)Nl~"Vr,{a) A i?f, 

N^oo 

where (-ff ) is given by \28\) . 

It is also possible to provide a multidimensional counterpart of Theorem [3] in the case of a vector of scales 
{ai)i<i<t where £ e N*. 

Theorem 4. Let {Xt)t>Q — {R^)t>Q be a Rosenblatt process and VN{a) be defined by f^. Then for every 
Q > 1 and LI G (i, 1) it holds that 



:7^^^Jv(a*)) A {R^„...,R^„ 

'Ct2{H) } \<i<l N^oo ^ ' ' 



where R^^ are normalized Rosenblatt random variables for all i = 1, ■ ■ ■ ,£ and for all Ai, . . . , G M the k-th 
cumulant of the random variable X]^=i ^j^i j 

ji,...,jk=i 

where Ck{R^) denotes the k-th cumulant of the Rosenblatt random variable with self- similarity order LL . 



are random variables in the 



Remark 1. Notice that, since the components of the vector (^R^^, . . . ,R^ f 

second Wiener case, its finite dimensional distributions are completely determined by the cumulants. 
Moreover, we deduce that the asymptotic covariance matrix of the random vector (N^^^VN{ai)] is 



i<i<e 



^,^CUH){{^3r"),^,,^r (29) 

Proof: The proof follows the lines of the proof of Theorem [31 Since for every Ai, . . . , Am £ R the linear 
combination Ui :— ^jTjfljT) ^Jv(aj) is a multiple integral of order two, it is possible to compute its 

cumulants by using the formula (I?!]) . We will obtain the following expression for the fc-th cumulant of Ui, 



Cfc(?7f)=( ^^°^^^ ) j^^ dyi ■ ..dyk{Xih^^\yi,y2) + ■■■ + ^r^^w"^ (2/1,^2)) 



X ... X (Ai/i^^(2/fc,2/i) + --- + 7j^^^"^(2/fe,yi)) 



15 



and then, as in the proof of Proposition [21 we can write 

/ c^wi / du[... duk duUT\\ui-u'i\\u'i-ui+i 

Jo Jo Jo Jo ^ 



where we used again the convention Uk+i = ui . We proceed as previously in the proof of Proposition [2] (with 
changes of variable) and then we obtain 



.2H'~2 



/ duidu'i . ..dukdu'J TT |M;a;; - - M;+ia;;+ij;+i + - ii+iii+i 



'[0,1]^'= j^i 

Thus, when Na — > oo, 



I G„fi/)C..fi/) ) E _7— — PtE---E j 



CkiUi) 



ji,...j"fc=l ^■^ ' ii=l ik = l 1 = 1 

e „ k 



— ^ 2J . . . Ajj^ Cfc(i?f ), 
where is a Rosenblatt random variable with index H and Ck{Ri) is its fc-th cumulant. 



Remark 2. is possible and instructive to study the behavior of the term T4 in the cases Q > 1 and 
H G (i, 1) or Q = 1 and H G (i, |). /t can be already seen from its asymptotic variance that it is very close 
to the Gaussian case. We can actually show this term converges in law to a Gaussian random variable. This 
fact does not influence the limit of the statistic Vn but we find that it is interesting from a theoretical point 
of view. We will denote by Ct^ (H) a positive constant such that 

NaETi -^N^^oo C^SH). 

Then the following holds: suppose that is a Rosenblatt process with self- similarity order H . Suppose that 
Q > 1 or Q — \ and H (£ (5, |)- Then with T4 defined in 



Nan A N{0,ClSH))- 



N 



The proof can be done by using a criteria in 124] in terms of the Malliavin derivatives. In the other case, 
i.e. Q = 1 and if G (|, 1), the limit in law of the renormalized Vm is a Rosenblatt random variable. 

5. Applications and simulations 

Throughout this section = (X/^)tgR^ denotes a fBm with H E (0, 1) or a Rosenblatt process with 

ffe (1/2,1). 
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5.1. Asymptotic normality of the sample variance oj approximated wavelet coefficients 

Here a sample [X^ , ■ ■ ■ ,X^) of X^ is supposed to be observed. For any couple (a, 6), define the 
following approximations of the wavelet coefficients d{a, b) and of the normalized wavelet coefficients d(a, b) 
defined in ^ and g]): 

e(a, b)^^y - b) and g(a, b) = ^^^^^^ , (30) 



The above expression are the usual Riemann approximations. Define also for a > 0, 



1 

VN{a} = — Y.{i'ia,^)-l). (31) 



Remark 3. These approximations of wavelet coefficients and their sample variance can be directly com- 
puted from data for all mother wavelet ip. In the particular case of a multiresolution analysis with or- 
thogonal discrete wavelet transform (that means a = '2?), the very fast Mallat's algorithm (similar to FFT 
for the Fourier's transform) can also be applied to obtain a different approximation 6^/(2-', /c) of d{2^,k). 
However, when j — )> cxi, e(2^,k) ~ eM(2"',fc). Indeed, let (j) the scaling function of a multiresolution anal- 
ysis satisfying J (j){t)dt = 1 and J^p^(t)dt — 1. Then it can be established that (see for instance ^2^1) 
eM{'2',k) = 2-^/2 / dtij{2-H - k) Y^^^i 4>{t - s)Xs = 2-^/2 X;fLi Xs j] dti/;(2-J(t + s) - k)(j){t). When 
j — >■ oo, since (j) is compactly supported on I, we obtain 'ip{2^\t + s) — fc) ^ ■0(2^^3 — /c) for s G /. Then 
eA/(2^,fc) ^ 2^^/2 ^^^(2^Js — fc) Jj (l)(^t)dt = e{2^ , k) because J (l){t)dt = 1; more precisely, a first order 

expansion gives 6^(2-', k) ~ e{2^ , k) + Op{2~^), with an approximation error which is negligible with respect 
to the approximation error computed in the following Lemma\^ The very low time consuming of the Mallat's 
algorithm together with a straightforward computation of e{2^,k) without multiresolution analysis, provides 
a clear advantage of the wavelet based estimator of the parameter H with respect to the estimators based on 
a minimization of a criterion (such as maximum likelihood estimators). 

Now the following result can be proved ; 

Lemma 1. Assume that ^ e C'"(M) with m > 1 and ip is [O,l]-supported. Let {ak)keN be a sequence of 
integer numbers satisfying N aZ^ — > oo and ajv — > oo. Then, for any Q > 1 , 

N^oo N^oo 

E\VMia^) - V^ia^)] < + + ^^_Qy,+H h2H-Q>-i)) ■ (32) 

^V"jv Oat a]^ " ^ 

Proof: First note that, 

1 / /"V^ 9 /-i ^ i.' / k' 



E[{iia,^)-d{a,^)r]=———( / / + - dtm H') (\t 



-t\'"-\t- 



2H 



a — 1 J 
k.k'=Q 



1. k ,2H 

U + -\ 4 
' a ' 


- \i H 


k' .2H 


|fc 


k' ,2H\ 




a ' 


' a 


a ' / 



^From standard Taylor expansion, if g is supposed to be m times continuously differentiable and [0, 1]- 
supported, for all a > 0, 



,^E^(^)-/ 9mt 



< sup \g^^^\t)\ 
te[o,il °' 



jl ij{t){i+tf" dt <C{{l+if"-'^). Therefore, 



and there exists C depending only on H , Q and tp such that 
there exists C depending only on H, Q, m and -0 such that for all a > 



C , |l^ , .fcM. k,2H 



I a ^-^ a a™ I a ^-^ a ' a' a™ 

k=0 ' " 



a'" I a ^ — ' a a 

k=0 
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Finally, as it was already proved in [3|, there exists C depending only on H and ijj such that for all m > 1 
and a > 0, 



a — ' a a a a 

k,k'=0 



,n2H 



C 

< — 

a 



< — . 

a 



fc'=0 

All those inequalities imply that there exists C depending only on H , Q, m and ip such that for all a > 0, 

E[(e(a,z)-d(a,z)) J <C[- + — + ^2,n j ■ (^3) 

Using Cauchy-Schwarz's inequality, 

N. 

n' 



^VN{a)-VN{a)\ < -^J2E\i^{a,i)~(F{a,i)\ 

< ^f]v'E[(e(a,*)-J(a,*))VE[(e(«,*)+rf(«,0)'] 

i—1 

1 ^" / / 

2 — 1 

1 1/2 1 ^° 

i—1 i— 1 

1 1/2 1 1/2 

^ ^(]^ E^[(^>'*)-'^>'*))']) (4 + 2— ^E[(e~(a,z)-d(a,z))']) 

i—1 i—1 

since E[(i'^(a,z)] = 1 by definition. Finally, from inequality (j33D and with a large enough, 

2—1 

/I log(iV) N^H~Q j^2H 

■ 

We will use Lemma [1] to prove the following result. 

Proposition 3. Assume that ijj G C™(M) with m > 1 and ^ is [0, l]-supported. Let {ak)ke'!Si be a sequence 
of integer numbers satisfying N aj} — >■ oo and aN — > oo. Then, 

N-^oo JV-i-oo 

1. ifX"isafBmwithQ<H<landQ>2orwithQ<H<3/4:andQ=l,andtfNa]^^~''^^^^ — > 0, 

then Theorem]^ holds when Vat (a) is replaced by Vat (a). 

2. if X" is a fBm with 3/4 < H < 1 and Q = 1, and if N a'^^^''^^^^ then Theorem holds 



when Vat (a) is replaced by Vat (a). 

f, 1 /O ^ ZT ^ 1 /O ^ 1 nf AT r, 

N 



3. if is a Rosenblatt process with 1/2<H<1,Q>1 and if N a^^ — > 0, then Theorems\^ and 



N^oa 

[2] hold when Vn (a) is replaced by Vn (a) • 

Proof of Proposition [3j The three different cases are respectively obtained from the Markov Inequality. 
Indeed, for e > and with a — 1/2, 2 — 2H, 1 — H (respectively), 

^ 1 A^" 

P{N:\VN{a)-VNia)\ > e) < - — E\Vn {a) - Vn (a)] . 
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Using Lemma [U it remains to obtain conditions for insuring 



a" \^/a^r a"^'" '-^+H / N^oo 



Uh-q>^i] (34) 



for all H to show the three cases of Proposition [31 

1. For a = 1/2, Q > 2 and < H < 1, condition ^ with a = 1/2 leads to 

max(iVa^2 ^ ^rQ^i-^ftFT ^ n a^f'^^""^ 12H-q>-i) < max(iVa^2 , A^a^'"^) 

It induces the condition N Uj^^ ^ ■* — > 0. 

For Q = 1 and any < H < 3/4, condition ([34]) leads to 

max(7Va^2^ A^a"^"^ , iVa^^"^) — ^ 0. 

N—i-cG 

Since < H < 3/4, it also induces the condition N a^^ ^ ' — > 0. 

2. if Q = 1 and 3/4 < < 1, then condition ([Mj) with a ^ 2 ~ 2H leads to 

_ 1 _ 1 1 _ 2m 1 m 1 m 

™^ C AT 4-4H ]\T ^ 4-2H ]\t ^ 3-2H -i AT" 2-2H -i \ 

max(A/ , Na^ , Naj^ Iq^i , iV lQ-2j 

< max(A^a^^ , N a^^~^) — ^ 0, 

and it also induces the condition N a,^^ 3 ) — ^ q 

JV-i-oo 

3. if Q > 1 and 1/2 < H < I, then condition ([34|) with a = l- H leads to 

max{Na'^''™ , iV a-i"™ , iV a^i-™lQ=i) < max (iV a^^' , iV a^'"") 



TV- 



which leads to the condition A^a^^ — > since to > 1. ■ 

Remark 4. Jbr a concrete estimation of H , the conditions between N and provided in Proposition\3[ 
do not depend on H . Usually, when this conditions depend on H , the convergence rate in the model can he 
improved. An adaptive procedure for estimating the smallest order possible for (ofc) could be also built as 
in the paper 0/. Anyway, we do not think that conditions provided on {ak)k in Proposition\^ are optimal. 
They could be improved by controlling E[(y/v(a) ~ VAr(a))^] instead o/E|V/v(a) — V;v(i)|- However, such 
computations are very long and tedious in the case of the Rosenblatt process ( it requires the computations of 
fourth-order moments) and we have preferred to avoid them. Moreover, the simulations we give below will 
show that our results are not so far to be optimal. 



Since the case of fBm has been already studied (see for instance [3|) we only provide below the numerical 
results when is a Rosenblatt process. Thus, we first exhibit the main result of this paper, i.e. the limit 
theorem C^^ {H) { ^) V/v(aAr) — > Ri following the procedure described in the sequel. 



Concrete procedure of simulations: 

• The samples of Rosenblatt processes are obtained following a similar procedure as the one presented 
in 0. With more details, to generate a Rosenblatt process sample {Xj^)i<j<N' 

1. generate a sample of length 1 + * to (in practice we use to — 100) of a fBm with parameter 
{H + l)/2. We use a wavelet based method introduced by Sellan (see [20j) with a Daubechies 
wavelet of order 10 but a circular matrix embedding method can also be applied (more details 
can be found in [3|)- Note that this sample is normalized and thus Var [/i?TO(l)] — 1. Next, one 
obtains a sample of length A^ * to of a fractional Gaussian noise (fGn) defined by the increments 
of the fBm, i.e. fGn{k) ^ fBm{k + l) - fBm{k) with Var/Gn(fc) = 1. 
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H 


0.6 


0.7 


0.8 


0.9 




[iVU.4] 






[AT"-*] 






[Aru.4] 






[Aru.4] 






N = 500 




3.31 


40.0 


3.35 


1.69 


40.8 


1.70 


1.08 


35.3 


1.49 


1.53 


114 


1.05 




IpMH 


2.23 


2.40 


2.11 


1.51 


1.55 


1.80 


0.85 


0.91 


0.94 


0.70 


0.75 


0.89 




i'c 


2.29 


2.45 


2.09 


1.76 


1.57 


1.49 


1.23 


0.85 


0.93 


1.05 


0.74 


0.81 


N = 2000 


V'4 


6772 


3.19 


2.99 


27935 


7.2 


1.64 


37951 


8.8 


1.56 


43888 


9.1 


0.75 




i'MH 


2.09 


1.91 


1.99 


1.33 


1.60 


1.40 


1.04 


1.08 


1.30 


0.68 


0.67 


0.77 




V'C 


1.95 


1.81 


1.72 


1.18 


1.41 


1.56 


1.21 


1.01 


1.25 


0.66 


0.72 


0.89 


N = 10000 




2.60 


2.77 


2.86 


1.49 


1.41 


1.43 


0.67 


0.84 


0.93 


0.52 


0.46 


0.59 




tpMH 


1.86 


1.79 


2.07 


1.10 


1.24 


1.43 


0.72 


0.70 


0.83 


0.51 


0.54 


0.60 




V'C 


1.55 


1.86 


1.72 


1.01 


1.08 


1.29 


0.64 


0.73 


0.75 


0.51 


0.53 


0.56 



Table 1: Computation for different choices of tp, H, N and ajv of y' MSE ofCj,^{H) ^Vjv(a]v) from 100 independent 

replications. 



2. Taqqu proved in [291 that: {^EtJlUGn^t) - 1)) A (ii-f )o<t<i, where i?^ is a (non- 

— — n— >oo 

normahzed) Rosenblatt process, since (1 + H)/2 > 3/4 is the parameter of the fGn. Thus, with 
n = mN, t = j/N and j ~ 1, . . . , N, one computes Yj — ^^j^^^ X]i='i(/^"^(*) ^ !)■ this way 
{Yj)i<j<N approximatively provides a path of (i?^^, Ryjf, ■ ■ ■ i Ri)- 

3. compute now = \J~^ff^^ Yj- for j = 1, . . . , A^: using the ff-self similarity of the Rosen- 
blatt process, (X/^, . . . ,X^) is approximatively a trajectory of a normalized {i.e. Va.r B.^ = 1) 
Rosenblatt process. 

The Matlab procedures to generate a trajectory of a fBm or a Rosenblatt process can be downloaded 
from |http: //seimos -univ-parisl -f r/-JeaxL-MeLrc-Bardet[ 

• Several mother wavelets "0 sue used: 

— The Daubechies' wavelet of order 4, (which is such that Q — 4 and C^([0, 1]) but not C^([0, 1])); 

— The Mexican Hat wavelet, ipMH (which is such that Q ~ 2 and E C°^(K) and is essentially 
compactly supported); 

— The function ipc such that Vc(i) = t{t - l){2t - l){t'^ - t + ^) for all t e [0, 1] and ipc = 
elsewhere (which is such that Q = 3 and tp e C°°{[0, 1]) except in and 1). 

• The values of constants C^{H), C,p{H') and Ct2{H) are obtained from usual approximations of inte- 
grals by Riemann sums. 

Montecarlo experiments using 100 independent replications of trajectories are realized for each H — 0.6, 0.7, 0.1 
and 0.9 and for N = 500, N = 2000 and N = 10000. The sequence of scales {aN)N is selected to be such 
that QN — [N'^''^], cln — [N'^'^] ot gn — [N'^'^]. The following Table [1] provides the results of simulations, 
which are values \/MSE of Cj^^{H) ( ^] V/v(aAf)) for the different choices of ip: H, N and gn- 

The main points to highlight from these simulations are the following: 

• globally, for TV and gn large enough then {C^^{H) ( ^ ) ^w(flw))jY seems to converge in distribu- 
tion to a centered distribution with a variance close to 1; 

• if Oat is not large enough (in these simulations, gm — [iV*^ "*]), there is a bias which clearly appears 
in the case of "04 since this wavelet function is not regular enough. In this way we may see that the 
conditions required for {gn) in Proposition [3] are close to be optimal. 

• Since Rosenblatt processes are generated from an approximation algorithm based on a wavelet synthe- 
sis, the larger H the higher the octave required to obtain a convenient trajectory. Because memories 
of computers are finite and we expect a low consuming time, the generator of Rosenblatt process has 
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Figure 1: FFT estimation (Silverman's method) of the density of the limit of {Cj,^{H) {^~~^ (o-n)) for H = 0.7, 

A'' = 10000 and fljv = A^"'^ from 100 independent replications. 




Figure 2: Convergence of the sequence {Cj,^(H) y-^'^ (.^n)) for H = 0.7 and 

a lightly smaller variance than it should have in the case H = 0.9. Moreover in such a case, there is a 
very slow convergence rate, i.e. (7V/aAr)°-^, in the limit theorem. 

An example of the estimation of the limit density is also presented in Figure[T]in the case H = 0.7, N = 10000 
and ttN = N'^'^. Such a density is quite similar to a standard Gaussian density but a Kolmogorov-Smirnov 
test invalides the hypothesis that this distribution is a Af{0^ 1) law. This result should be compared with the 
numerical simulation of the Rosenblatt density given in }3l| . 

Finally, Figure [2] shows the convergence of the sequence (C^^ (H) i-^j ^AfC^Af))^ when TV increases. 
This sequence does not seem to converge in L^(f7) as it is claimed in Section [3] 



5. 2. Estimation of H 

Here we consider that a sample {Xi, . . . , Xn) oi X = X^ with X^ a fBm or a Rosenblatt process is 
known and H and cr^ > are unknown. Denote the sample variance of wavelet coefhcients 

1 

Then, following Proposition |31 one deduces that 
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where the asymptotic distribution (£i)i<i<£ is a Gaussian distribution (with a = 1/2 when is a fBm 
and Q > 2) or a Rosenblatt distribution as defined in Theorem 4 (with a = 1 — H when is a Rosenblatt 
process). Therefore, from the so-called Delta- Method, we also have 

N:(log{TNiiaN))-i2H + l)\og{ia^)~\og{a^C^{H))) A (£.)i<.<^- (36) 

\ / l<i<£ W— >cxD 

Therefore, a log-log-regression of (/Ar(iaAr)^ by (*aA')]^<j<^ provides an estimator of H (this estimation 
method has been introduced in [is'l). Such an estimator is defined by 

Hn ■■= {\ , 0)' • {Z',Zey'z',{\og{lN{taN))),^^^, - \, (37) 
where Zi{i^ 1) = logi and Zf(i, 2) = 1 for alH = 1, • • • Then Proposition |3] implies 

Proposition 4. Assume that -0 G C'"(M) wit/i to > 1 and is [0, 1]- supported. Let {ak)ke'N be a sequence 
of integer numbers satisfying N ajj^ — > oo and — > oo. Let (Xi, . . . , Xn) an observed sample of 



X ^ X^ where X^ is a fBm or a Rosenblatt process. Then, 
fB 

then with (H , i/j , £) > defined in |^Q[ ) and depending only on H, xp and i, 



l.ifX"isafBmwithO<H<landQ>2orwithQ<H<3/4:andQ = l,andtfNaj^^ ^^^'"^ 0, 



— {Hm-H) a AA(0,7'(i?,V',^)); (38) 
2. if X^ is a Rosenblatt process with 1/2<H<1,Q>1 and if N — > 0, then 



N- 

{—Y~"{H^-H) A iH,.^., (39) 

where Ln.ip.e defined in |^_?[ ) is a distribution depending only on H , -0 and H.. 

Proof of Proposition [4) From the results of Proposition [3J the relation ([55)1 is clear since In {a) = 
Vat (a) + 1 for all a > 0. Now using the usual multidimensional Delta-method (see for instance 3J|) with 
the transformation function g{xi, . . . ,xg) = ( log(xi), . . . , log(a;f ))' applied to the limit theorem psp . one 
obtains _ 

where Jg(l, . . . , 1) is the Jacobian matrix of g at point (1, . . . , 1). Therefore, since Jg(l, . . . , 1) is the identity 
matrix, one obtains (1361) . Then, 2H + 1 can be estimated from an ordinary least square regression and we 
obtain: 

2Hn + 1 = (1,0)' ■ (Z;z,)-'z; (log(/jv(za^)))^<^<^ - Af, (log(/Ar(zajv))) 
with Me a (1 X ^) matrix. Hence, the formula ([57]) can be deduced. Moreover, (1551) implies 



n: 



r„ ((2 Hn + 1)- {2H + 1)) A M, (£.)i<K^- 



l<i<e N^oD 

Therefore, when X^ is a fBm, since (ei)i<i<f ~ 7V(0, Lf\H)) from dH]), we deduce (1381) with 

7^(i?,V,^) :=^M,Lf)(iJ)M;. (40) 
The same trick implies ([39)) when X^ is a Rosenblatt process and 

iff,^,£:=^M,(i?fi,...,i?f^,)' (41) 

with the distribution of [R^i, ■ . ■ , Rif) defined in Theorem |31 H 
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Remark 5. An estimator of can also be provided by this method, with the same convergence rate. 

Remark 6. In Proposition^ we do not study the case where is a fBm, Q = 1 and 3/4 < H < 1. 
Indeed, in a statistical framework devoted to the estimation of H , since the choice of tp is arbitrary, there is 
no reason to chose tjj with Q — 1 which gives a worst convergence rate than tp with Q > 2. 

We summarize our results concerning the convergence rate, with m > 2 and ai\i ~ ]\[^/'^+^ with S > 
arbitrary small: 

1. if is a £Bm and Q > 2 or Q = 1 and < < 3/4, the convergence rate of Hn is 7Vi/4-<5/2. 

2. if is a Rosenblatt process and Q > 1, the convergence rate of Hn is A''^-'^"^^/^"*^^"^^ . 

Such convergence rates are weak in a parametric framework. For instance, applied to the increments of 
a fBm, the convergence rate of the maximum likelihood or the approximated Whittle maximum likelihood 
estimator are N^^'^ (see 11| or [3]); as far as we know, there are not such results in the case of the Rosenblatt 



process. An estimator based on quadratic variations has been studied in [33] and a noncentral limit theorem 
is proved with convergence rate N^^^ ; but such an estimator is almost parametric and can not be applied 
to processes which are not strictly self-similar. As it was previously recalled in the introduction, the wavelet 
based estimator is interesting because it can also be applied in a semiparametric framework. 

Remark 7. In a statistical framework when X^ is a Rosenblatt process, {Xi, ■ ■ ■ ,Xn) is known but H 
unknown, how to obtain confidence interval from ^39\) ? We propose a three step procedure. The first step is 
to chose On = [iV^/^"*""^] with 6 > arbitrary small (for instance S — 0.05) and compute Hn- The second step 
consists in the computation of the quantile q (for instance the 97 .5%i-quantile) of the distribution of L^^ ^ ^ 

(defined by i39\)) from Monte-Carlo simulations of Rosenblatt processes with parameter H = H^. Indeed, it 
is possible to replace I139\} by: 

{Hn - H) - L^^^^ i, 0, 



because L^^^^^^ - Lh^^^i and {^f x {^)" ' = exp [{H - HN)\og{^)] 1 from 

h39\) and Slutski's Lemma. Finally, the confidence interval of H will 6e [ — g (^)^" 9 ^] since 

the density of L ^ i is clearly an even function. 

Since there are not other asymptotic results than ours on semiparametric estimators of H for a Rosenblatt 
process, we can only numerically compare the wavelet estimator with other estimators. We have chosen two 
usual semiparametric estimators of long memory parameters: the local Whittle estimator (denoted Hwhutie) 
defined in ^26J and the adaptive global log-periodogram estimator (denoted HLogPer) defined in [2^ . Remark 
that such estimators are applied to the increments of X. Note also that both the Matlbab softwares of these 



estimation procedures can be downloaded from http://samm.univ-pcirisl.fr/-Jeaii-Marc-Bardet For 



completing the numerical study, we have considered two cases of wavelet based estimators: 

• Hnaar whcrc ip is the Haar mother wavelet defined by ip{x) = 1 for < x < 1/2 and iIj{x) ~ —1 for 
1/2 < a; < 1 and elsewhere ip — 0; moreover, we consider un — [N'^'^] and £ — [N'^-^]. 

• where ipc is wavelet function defined above; moreover, we consider un — [A^°'^] and £ — [N'^-^]. 

For both these estimators the number £ of scales is chosen as a sequence depending on N because additional 
simulations have shown that the variance of these estimators is optimal with such a choice. The results of 
simulations are given in Table [521 

The conclusions from Table 15.21 are the following: 

• The local Whittle estimator Hwhittie is clearly the most accurate. The wavelet based estimator H^^ is 
almost as accurate as Hwhittie- The estimator HLogPer seems to be a little less accurate while Hnaar 
is not satisfying (it can be explained by the fact that this wavelet function is not continuous in 1/2). 
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H 


0.6 


0.7 


0.8 


0.9 






Mean 


VMSE 


Mean 


VMSE 


Mean 


VMSE 


Mean 


VMSE 


N = 500 


— 


0.38 


0.32 


0.43 


0.35 


0.47 


0.37 


0.43 


0.49 




Hi, 


0.56 


0.15 


0.62 


0.16 


0.72 


0.15 


0.84 


0.10 




W ritttLe 


0.62 


0.07 


0.69 


0.07 


0.77 


0.08 


0.87 


0.09 




if T ^ ^ D ^ « 

^-^ LlOg r er 


0.58 


0.17 


0.67 


0.16 


0.71 


0.18 


0.77 


0.24 


N = 2000 


^ H aar 


0.50 


0.18 


0.60 


0.18 


0.70 


0.16 


0.83 


0.11 






0.59 


0.07 


0.65 


0.07 


0.75 


0.08 


0.85 


0.07 




^Whittle 


0.65 


0.06 


0.71 


0.04 


0.80 


0.05 


0.88 


0.05 




H LogPer 


0.60 


0.08 


0.65 


0.10 


0.75 


0.12 


0.82 


0.13 


N = 10000 


^Haar 


0.54 


0.16 


0.58 


0.15 


0.66 


0.16 


0.69 


0.22 






0.60 


0.04 


0.66 


0.05 


0.76 


0.05 


0.86 


0.04 




^Whittle 


0.64 


0.04 


0.72 


0.03 


0.79 


0.03 


0.89 


0.03 




^ LogPer 


0.60 


0.05 


0.67 


0.06 


0.74 


0.08 


0.85 


0.07 



Table 2: Empirical mean and V MSE of Hnaar, H^^, H^f^ittle, HLogPer for different choices of H and N, from 100 
independent replications. 



• It appears that the convergence rate of the estimators depends on N and all the estimators seem to 
be consistant. 

• The variance of HLogPer seems to increase as N^~^ when H increases unlike both wavelet estimators 
while theoretical results say that their variances behave as N^~^ (see Proposition |4]). An explanation 
of this phenomenon is the following: the asymptotic variance Var(iJ^) of '^^ Hnaar behaves as 

the (2, 2)-component of the matrix — (Z^ Zij Z'^ Zg (Z^ Zg^ ( ^ j where is given in 
(P^. After computations we obtain that 

Therefore the larger £ the smaller Var(if^). Moreover if TV = 10000, ajv = [A^°"] and £ = [N°-^], 
with ?/> = ipc, we obtain for Va.T {%) ~ (0.024)2, (0.028)^, (0.030)2 ^^^^ ^ (0.030)^ for respectively 
H ~ 0.6, 0.7, 0.8 and 0.9. It explains why the different vMSE do not numerically seem to depend on 
H even if they theoretically depends on H. 
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